International Journal of Engineering Research & Science (IJOER) 


ISSN, [2395-6992] 


[Vol-2, Issue-5 May- 2016] 


Stability and Hopf bifurcation for Kaldor-Kalecki model of 
business cycles with two time delays 

Xiao-hong Wang 1 , Yan-hui Zhai 2 , Ka Long 3 

School of Science, Tianjin Polytechnic University, Tianjin 

Abstract — Papers investigate a Kaldor-Kalecki model of business cycle system with two different delays, which described 
the interaction of the gross product Y and the capital product K . We derived the conditions for the local stability and the 
existence of Hopf bifurcation at the equilibrium of the system. By applying the normal form theory and center manifold 
theory, some explicit formulate for determining the stability and the direction of the Hopf bifurcation periodic solutions are 
obtained. Some numerical simulations by using Mathematica software supported the theoretical results. Finally, mam 
conclusions are given. 
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I. Introduction 

According to rational expectation hypothesis, the government will take into account the future capital stock in the process of 
investment decision. Business cycle, named economic cycle, refers to the total alternate expansion and contraction in 
economic activities, and these cyclical changes will appear in the form of fluctuation of the comprehensive economic activity 
indicators such as gross national product, industrial production index, employment and income. In the early research of non- 
linear business cycle theories, Kaldor [1] assumed that investment depended on gross product and capital stock and proved 
that, through graph analysis, time-varying nonlinear investment and saving function could result in a business cycle. Chang 
and Smyth [2] summarized Kaldor's theory on business cycle, established a nonlinear dynamic system which described the 
gross national product and capital stock changed over time, gave necessary conditions of existence of limit cycles, and 
provided a rigorous mathematical proof of Kaldor's model proposed in 1940. On this basis, Grasman and Wentzel [3] further 
improved Kaldor's business cycle model by considering capital loss speed. Along another view, according to the IS-LM 
model raised by J.R. Hicks and A.H. Hansen, Ackley [4] established a complete Keynes system which reflects the gross 
domestic product and interest rate changes over time, which is also called standard IS-LM model. Kalecki [5] first considered 
investment delay in the business cycle model in 1935, and he claimed capital equipment needed a conceived cycle or delay 
from installation to production. As the theory of delay functional differential equations gradually become more accomplished 
in 1990s, Krawiec and Szydlowski [6,7] first made a qualitative analysis of the impact of the investment delay on the 
business cycle. In addition, some other scholars investigated the dynamic properties like stability, many types of bifurcation, 
existence and stability of periodic solutions in Kaldor-Kalecki model with investment delay [8,9]. The record of business 
cycle has been kept relatively well during the last 200 years, and business cycle theory, as the core issue of macroeconomics, 
has been attracting the widespread interests of many economists. The modern business cycle theory can be traced back to the 
masterpiece of Keynes's theory "The General Theory of Employment, Interest, and Money". Keynes discussed the formation 
of the business cycle from the perspective of psychological factors based on national income theory. The following system 
was formulated by Krawiec and Szydlowski [10] who combined two basic models of business cycle: the Kaldor model and 
the Kalecki model. Kaldor [11] first treated the investment function as a nonlinear (s-shaped) function on Y so that the 
system may create limit cycles, while Kalecki [12] assumed that the saving part is invested, and that there is a time delay due 
to the past investment decision. Therefore, the gross product is available in the market after a time lag. 

Y(t) = a(I (Y , K) - S (Y , K)), 

< 

K(t) = I(Y(t-r),K)-SK. (11) 

Y is the gross product and K is the capital product of the business cycle; (X > 0 measures the reaction of the system to the 
difference between investment and saving; 8 e (0,1) is the depreciation rate of capital stock; I, S '. R x R — > R are 

investment and saving functions of Y and K , respectively; T is a time lag representing the delay for the investment due to 
the past investment decision. In a business cycle system, delay occurred in the production of investment function not only, 
also released on capital stock, so to introduce the following model. 




International Journal of Engineering Research & Science (IJOER) 


ISSN, [2395-6992] 


[Vol-2, Issue-5 May- 2016] 


Y(t) = a(I(Y , K) - S (Y , K)), 

K(t) = I(Y(t-T),K(t-T))-SK. (12) 

However, for the systems (1.2), the delay that occurs in investment function and happened in the capital stock is not 
necessarily the same, so we set up two different time delay. By quoting the function of investment and capital saving 
function in reference [1], I (F, K ) = I (F) — (3K , S (F, K) — yY , f3 > 0 , y e (0,1) . Following model be produced. 

Y(t) = a(I(Y)-j3K-yY), 

< 

k(t) = I(Y(t -T x ))-pK(t-T 2 ))-SK. (1 3) 

The rest of this paper is organized as follows. In Section 2, we analyze the distribution of eigenvalues of the linearized 
system of (1.3) and derived the conditions for the local stability and the existence of Hopf bifurcation at the equilibrium of 
the system. In Section 3, by applying the normal form theory and center manifold theory, some explicit formulate for 
determining the stability and the direction of the Hopf bifurcation periodic solutions are obtained. In Section 4, some 
numerical simulations by using Mathematica software supported the theoretical results. Finally, main conclusions are given. 
In Section 5, main conclusions are given. 

II. Stability and Analysis of Local Hopf Bifurcation 

Let E*—(Y*,K *) be an equilibrium point of Sys.(1.3), I*—I(Y *), and u l =Y—Y*, u 2 — K — K *, 
/(.S') = I(s + F*) — I* . Then Sys.(1.3) can be transformed as 


u x = a(i(u x )- yu l — pu 2 ), 
u 2 = i(u l ( t - r, ))- Pu 2 ( t - t 2 )) - Su 2 . 

Let the Taylor expansion [13] of / at 0 be 

/(«) = I (F*) + ^ - (u) 2 +o(\u I 3 ). 

Then Sys.(2.1) can be transformed as 

• , t /"(f*) 

= a{l (F*)m, 4 — (wj) 2 -yu x - pu 2 ), 

• I (F*) 

u 2 = I (Y*)u x (t -T x ) 4 (uPt-Tj) 2 - Pu 2 (t-T 2 ))-8u 2 . 

Then the linear part of Sys.(2.2) at (0,0) becomes 

u x = a(I (F*)mj -yu x - Pu 2 ), 
u 2 =/ (Y*)u x (t -Tj)- Pu 2 (t -t 2 ))-8u 2 . 

And the corresponding characteristic equation is 

A 2 4- (8- - a) X + (A - a) pe^ T2 + be^' -8a = 0, 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


Page | 2 




International Journal of Engineering Research & Science (IJOER) 


ISSN, [2395-6992] 


[Vol-2, Issue-5 May- 2016] 


where a = a {I (F*) - y)A,b = a(5I (F*). 

Theorem 2.1 When r, = r 2 = 0 , the characteristic equation of system (2.4) is 

A~ + (r?> — ci + + b — §ci — Pci = 0. (2.5) 

All roots of the equation (2.5)has negative part if and only if 

(H n )A l0 = S-a + fi > 0, A u =b — 5a— /fa > 0. 

Then the equilibrium point E* — (Y * , K ) is locally asymptotically stable. 

Theorem 2.2 When r, > 0 ,T 2 = 0 , assume that (H n ) is satisfied. Then (2.2) has a pair of purely imaginary roots + i(0 U) 
when T = r 10 . 

Proof.The characteristic equation of system (2.4) is 

A 2 + (8 - a + p)k + be~ Az ' - a(S + p) = 0. (2.6) 


Clearly, id ) , is a root of Eq. (2. 6) if and only if i(0 { satisfies 

-cof +{S-a + P)co l i + b{co?,co x z x - i sin co^p-aiP + d) = 0. 

Separating the real and imaginary parts, it follows 

j- 0 ) 2 +bcosco l r 1 =a(P + S), 

[(^-a + P)co x =Z7sin<z) 1 r 1 , 


(2.7) 


According to conclusions of Beretta and Kuang [14], the stability of the system changes, when the real part of its 
characteristic root passes through zero point. Therefore, considering the critical situation, let characteristic root A real part 
Re A = 0, , adding up the squares of both Eq.(2.7), it yields 


vf +(c + a 2 )v l +a 2 c-b 2 = 0, (2.8) 

where Vj = CO 2 ,C = {S + /?)" . 

Hence Eq.(2.8) has solution V 1() and v n .Where 

-( a 2 +c) + J(a 2 +c) 2 -4 (a 2 c-b 2 ) — ( a 2 +c)~J(a 2 +c) 2 -4(<7 2 c-& 2 ) 

V|0_ 2 ’ V ‘ 1_ 2 

Assume 

h(V\ ) = vf +(c + a 2 )v 1 +a 2 c-b 2 =0, (2.9) 

Noting that v n is negative. Eq.(2.9) has one positive root if and only if (// 0| )!y > Cl~C . The characteristic Eq.(2.6) has a 
pair of purely imaginary roots + ico l0 = +i-Jv w . And, hence 




O, 


10 


The proof is completed. 
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Differentiating both sides of (2.6) with respect to T ] , it follows that 


dA , _ 2A + d + /5-a r, 
dr t -A be~ Zr ' A 


That is 


p , dA , | A=iC0l0 _ 2 co 4 + ((d + [5) + a~)ca~ 

^ ' * i = r 10 / _/,\2 


d r, 


{coby 


Obviously, it is greater than zero. So, the transversality condition is meted. By applying Lemmas 2.1 and 2.2 and condition 
(#2l) , we have the following results. 


Theorem 2.3 For system (2.2), if the condition is satisfied, the equilibrium E of system (2.2) is asymptotically stable for 
Tj G [0, r 10 ) . System (2.2) exhibits the Hopf bifurcation at the equilibrium E for f, = r |0 . 


Theorem 2.4 When T x = 0 , Z" 2 > 0 , assume that {H M ) is satisfied. Then (2.2) has a pair of purely imaginary roots + ico 20 
when T = . 

Proof.The characteristic equation of system (2.4) is 

A" + ( 8 — ci) A + (A — ci) fde 1 + b — ciS = 0 . ( 2 . 10 ) 


Let id) <2 ( CO > 2 > 0) is a root of (2.10). Then 

-cv 2 +(S -a)co 2 i + j3(ico 2 cos co 2 t 2 +ia sin co 2 t 2 +co 2 sin cq 2 t 2 -acosoo 2 T 2 ) + [5 -ad = 0. 


Separating the real and imaginary parts, it follows 

j-co 2 + (dco 2 sin co 2 t 2 -a/3cosa> 2 T 2 -ad + b = 0, 
[(S - a)a> 2 + ( 3 < x > 2 cos co 2 t 2 + [5a sin a> 2 r 2 . 


( 2 . 11 ) 


Adding up the squares of both Eq.(2.11), it yields 

v 2 + (S 2 - P 1 + 2a 2 - 2 b)v; - (2 abd -2 a 2 (d 2 - J3 2 )- (a 2 - bf)v 2 + 
a 2 b(b -lad) + a 4 (d 2 - /3 2 ) = 0, 


where 0) 2 = V, . 

From the analysis above, we can obtain that if all parameter values of the system (2.2) are given, the root of the equation 
(2.12) can be easy to solved by use of the Mathematica software. Therefore, in order to give out the main conclusion in this 
section, the following assumptions is given. 

(// 31 ) The equation (2.12) at least has one negative root. 


If the condition (//,,) is satisfied. Then, the characteristic Eq.(2.12) has an positive root V 20 .So, the Eq.(2.10) has a pair of 
purely imaginary roots + i(O 20 = ±z\Zv-, 0 . And, hence 


r 20 = — arccos(- 


ox 


20 


ab-d(a n +ty 20 ) 
P(a~ + co 20 ) 


)• 


The proof is completed. 
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Differentiating both sides of (2.10) with respect to T-, , it follows that 


dA _ 2A + 8 + fie AT2 - a r 2 
dr 2 PA(A - a)e~ Az2 A. 


That is Re( ) 1 \ A 


Jr, 


l0J 2C 

T 7 =T ?n 


2(d\ 0 + (J 2 + p 2 +3a 2 )cc>2 0 + a 2 (S 2 + p 2 + a 2 )a>l 0 
(fJ 2 co 4 + (afJS) 2 )(a 2 +co 2 ) 


Obviously, it is greater than zero. So the transversality condition is meted. By applying Lemmas 2.4 and condition ( /Y, , ) , 
we have the following results. 

Theorem 2.5 For system (2.2), if the condition (J/ 31 ) is satisfied, the equilibrium E of system (2.2) is asymptotically 
stable for r 2 €E [0, r 20 ) . System (2.2) exhibits the Hopf bifurcation at the equilibrium E for r, = r 20 . 


Theorem 2.6 When Z", = r 2 = T , assume that ( H 4t ) is satisfied. Then (2.2) has a pair of purely imaginary roots + i(D Q 
when T = r 0 . 

Proof.The characteristic equation of system (2.4) is 

A 2 +(S- a) A + {{A-a)p + b)e Ar -aS = 0. (2.13) 

Let i(0{C0 >0) is a root of (2.13). Then 


- cc>l+(8-a + Pe /T )coi-(aP -b)e At -a§ = 0. 

Separating the real and imaginary parts, it follows 

\-a> 2 + P<x>sm<x>z ~{aP~ Z?)cosor -ad = 0, 

[(S - a )(0 + Pcocoscot + ( Pa - b) sin ® 2 r 2 = 0, 


(2.14) 


v 3 + ( J 2 - P 2 + a 2 + 4 - a) 2 )v 2 + (2(^ - 8){^ - a)aP - 2 {ap~bf + 


Adding up the squares of both Eq.(2.14), it yields 

,b_ 

P 

) + . 

P P 


(a(a-^) + S^) 2 )v + (a-^) 2 ((ap-b) 2 ~(aS) 2 ) = 0, 


(2.15) 


where CO —V. 

From the analysis above, we can obtain that if all parameter values of the system (2.2) are given, the root of the equation 
(2.15) can be easy to solved by use of the Mathematica software. Therefore, in order to give out the main conclusion in this 
section, the following assumptions is given. 

(// 41 ) The equation (2.15) at least has one negative root. 


If the condition (H 4l ) is satisfied. Then, the characteristic Eq.(2.15) has an positive root V 0 
purely imaginary roots + i(X) Q = And, hence 


So, Eq.(2.13) has a pair of 


Aarccos ^bSP) + aS{b-aP) ) 

co Q P a>l +(a,p -b ) 2 
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The proof is completed. 

Differentiating both sides of (2.13) with respect to T , it follows that 

<tt , _ 2/1 + 8 + pe Ar - a r 
dr A(j3(A-a) + b)e~ Xr 2 

That is 

Re A- 1 \ A=ia> « - 2^0 +(2 (b-a/3) 2 + /3 2 (a 2 + 8 2 ))co 4 0 +(b-a/3) 2 (8 2 +a 2 )co 2 
far T=To (Jd 2 col+(aP-b) 2 ) 2 

Obviously, it is greater than zero. The transversality condition is meted. By applying Theorem 2.5 and condition ( H 4] ) , we 
havethe following results. 

Theorem 2.7 For system (2.2), if the condition (H r ) is satisfied, the equilibrium E of system (2.2) is asymptotically 
stable for T e [0, r 0 ) . System (2.2) exhibits the Hopf bifurcation at the equilibrium E for T = r 0 . 

Theorem 2.8 When T l > 0, r 0 e (0,r 20 ) , assume that (// 51 ) is satisfied. Then (2.2) has a pair of purely imaginary roots 
+ ico x when r = r*. 

Proof. Let i COj ( (0 t >0) is a root of (2.4). Then 

-fayj) 2 + (S - a)0J t i + (d(ico x -«)(cos&\r 2 -isinco l T 2 ) + b(cosco l T l -ismcd^-aS = 0. 

Separating the real and imaginary parts, it follows 

— (<z>,) 2 + Jdco, sintaV, -aBcosoo,T 2 + Z?cos<ur, -aS = 0, 

i r z ii (2.16) 

[( 5 -a )co x + (dco x cosfo»,r 2 + /fa sin -bsina> l T l = 0, 

Adding up the squares of both Eq.(2.16), it yields 

{co x ) A -(2/?sin co 1 t 2 )(cl> 1 ) 3 +(J3 2 + 8 2 +a 2 + 2/?£cos&> 1 r 2 )(fa 1 ) 2 - co x 2a 2 J3 sinryjr 2 

+ a 2 (j3 2 + 8 2 ) + 2a 2 J3 8 cos co x t 2 -b 2 = 0. (2.17) 

From the analysis above, we can obtain that if all parameter values of the system (2.2) are given, the root of the equation 
(2.17) can be easy to solved by use of the Mathematica software. Therefore, in order to give out the main conclusion in this 
section, the following assumptions is given. 

( H 5l ) The equation (2.17) at least has one negative root. 


If the condition (i/ 51 ) is satisfied. So, the characteristic Eq.(2.17) has an positive root i 0)[ . Eq.(2.4) has a pair of purely 
■ • 

imaginary roots ± l CO x . And, hence 


* i 1 

r, = — arccos(- 


ty, 


(2 a-S){o\) 2 +a 2 8-abcosa\r 2 +ba\ sinty*r 2 
{co x ) -a 


). 


The proof is completed. 

Differentiating both sides of (2.4) with respect to T x , it follows that 
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-A x ? 


_ 2X + 8 + pe 2 -a-T 2 P(h-a)e r x 


w - 


Abe 


-At, 


A 


That is 


dA ^ | lX=io f _ 2co* cos a\T* + (J3 — a) sin co*r* + fi( \ + ar 2 )sinfc» ! '(r | : -t 2 )~ t 2 Pco* cos co* ( t* —t 2 ) 


Re (— ) 1 1 
dr j 


0){h 


dA 


Obviously, Re( ) I ^ O.The transversality condition is meted. By applying Lemmas 2.8 and condition (H ) , we 

</r, ri “ ri 

have the following results. 

Theorem 2.9 For system (2.2), if the condition r 2 e [0, Z" 20 ) and (H 5l ) is satisfied, the equilibrium E of system (2.2) is 
asymptotically stable for r l E [0, rj*) . System (2.2) exhibits the Hopf bifurcation at the equilibrium E for T x = T* . 

Theorem 2.10 When T 2 > 0, T) e [0, r 10 ) , assume that (H h] ) is satisfied. Then (2.2) has a pair of purely imaginary roots 


+ ico 2 when Z 2 =t 2 . 


Proof.Let + ico 2 (a)* >0) is a root of (2.4). Then 

-( co 2 ) 2 + (d -a)co 2 i + P(ico 2 - a)(cosco 2 r 2 -(sinty 2 r 2 ) + Z?(cos<y 2 r 1 —i sin co 2 T x )-a8 = 0. Separating the real 
and imaginary parts, it follows 

\-(co P 2 + fico, sin^r, -a/?costy 2 r, + Z?cos<z> 0 z\ -ad = 0, 

i , , , (2.18) 

[(S-a)co 2 + pco 2 cos co 2 t 2 + fda^ma) 2 x 2 -bsm<x> 2 X\ = 0, 

Adding up the squares of both Eq.(2.18), it yields 

( co 2 ) 6 +(2 a 2 - p 2 +d 2 - 2b 2 cosry 2 rj)(ry 2 ) 4 -(2b(a-S)sin(d ] T 2 )(a) 2 ) 3 + 

(2a 2 (S 2 ~p 2 ) + a 4 +b 2 - 2ab(S + a) cos co 2 tP(co 2 ) 2 +2a> 2 a 2 b(a -5)sma> 2 T l + (2.19) 

a 2 (S 2 - p 2 +b 2 - 2ab 5 cos co 2 t x ) = 0. 

In order to give out the main conclusion in this section, the following assumptions is given. 

(H (it ) The equation (2.19) at least has one negative root. 


If the condition (H (j] ) is satisfied. Then, the characteristic Eq.(2.19) has an positive root (0 2 .Eq.(2.4) has a pair of purely 
imaginary roots + ico 2 . And, hence 


* 1 

t 2 = — arccos(- 

(Q-, 


(2 a-d)(a\) 2 + a 2 8 -ab cos co* 2 t x + ba> 2 sina> 2 T 1 
(a> 2 ) -a 


)■ 


The proof is completed. 


Differentiating both sides of (2.4) with respect to r, , it follows that 

dA j _ 2 A + 8 + Pe lT ' + Pe lH -a-Tpe~ 2T ' r 2 
dr Ap(A-a)e^ 2 ^ A 


Page | 7 
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That is Re(— ) 1 \ A ~ io j 

dz , r2=r2 


P\co 2 ) 2 {{co’ 2 ) 2 +a 1 y 


, where 


A- Pzol{a{5 -a)-2{oP 2 y)sma>* 2 z*2 - P (ry 2 ) 2 (S + a) cos (o* 2 z* 2 + 
-T^ + asinry*^ ~ T \))~ PS 0 ^) 2 ■ 


Obviously, Re( ) 1 ^ () . So, the transversality condition is meted. By applying Theorem 2.10 and condition 


dr , 


r 2 =r 2 


(H 6l ) , we havethe following results. 


Theorem 2.11 For system (2.2), if the condition r, e [0, f l0 ) and (H 6l ) is satisfied, the equilibrium E* of system (2.2) is 
asymptotically stable for Z~, E [0, z) ) . System (2.2) exhibits the Hopf bifurcation at the equilibrium E* for Z 2 = Z*. When 
T~, > 0, Tj e[0,r 10 ), assume that (H (] ) is satisfied. Then (2.2) has a pair of purely imaginary roots +ico* when 
r 2 =r*. 


III. Direction and Stability of the Periodic Solutions 


In this section, by applying the normal form and center manifold theory [15,16], we discuss the direction and stability of the 
bifurcating periodic solutions for Z 2 > 0, 7) E [0, Z" l0 ) . Throughout this section, we always assume that system (1.3) meets 

the conditions of the Hopf bifurcation. In this section, we assume that 7" | *<Z\,Z" 1 e[0, Z" 10 ). For convenience, let 

r 0 = jU + Z*, clearly, // = 0 is the Hopf bifurcation value of system (1.3). Let t / Z and still denote t .Then the system 
(1.3) can be rewritten as 

u(t) = L M + F(ju,u t ), (3.1) 

where u(t ) = ( U l (t),U 2 (t)) T and 


PSP) = « + + B'<K- % + £</>(- 1)), 

r 2 


f 

F(/u,</>) = (r 2 + jj) 

v 



where 


% 

1 

a 

, B = 

2 0 

0^ 

, and C = 

o 

o 

10 -S) 


PI 

Oy 


1 

o 


And q = I (Y*),F l = aqzjp (0), F 2 = q (j) 2 ( Z p L ),q = 1 ( ^ - ■ 

z 2 2 

By Riesz's representation theorem, there exists bounded variation functions 

r/(0,/u ) : [—1,0] — > R x ~ for 9 e [—• r*,0) ,such that 

LJ) = | d rj{0, /d)<f>(Q), <t> e C. 


In fact, we take 


(3.2) 


(3.3) 


Page | 8 
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\t 2 + ju)(A +B + C), 

9 = 0, 


(t 2 + ju)(B +C), 

9e[-^, 0), 

-7- 



~2 

(3.4) 

(t 2 + p)C, 

0e(-l,-%), 



Tl(0,fi) = 


[ 0 , e=-\. 

For I// G C [—1,0] , we define that operator A and R is as follows form. 

d(m) 

A(ju)m= 


d0 


9 e [—1,0), 


\d(q(0,p)mi 9 = 0. 

[0, 9 g [—1,0), 

R(»)m = \ 

1 FQi,0), 9 = 0. 

Then the system (3.1) is equivalent to the following abstract differential equation 

u t = A(ju)u t +R(ju)u t , 

where u, = + 9),u 2 (t + 9),u 3 (t + 9)) for 0e[-r 2 °, 0]. 

Setting I// G C [0, l\ , we can define an adjoint operator A*(0) corresponding to A(0) as the following form: 

dy/(s) 


A*y/(s) = 


ds 


s g (0,1], 


f° 0 d(7 r (s,0)^(-s)),s = 0, 

J T') 


(3.5) 


(3.6) 


(3.7) 


For cp G C ([— ■ T 2 ,0], R ) and y/ e C [0, r 2 ] , define the bilinear form: 

< <p(s),</>(9) >= <p(0)</>(0) - [ f <p(g - 9)d q(9)^)d ^ , (3.8) 

J- 1 J£=0 

where ?](9) = J](0.0). After discussed above, we know that + T*(0* is characteristic root of A(0) and A* . To determine 
the normal form of operator A , we need to calculate the eigenvectors p(9) = (1, p 2 ) T e r2< ° 2<> and 
p* (9) = D(l, p*) T e' Tl< ° 2S of A and A* corresponding to + T*C0* and — V 2 C0 * , respectively. By the definition of the 
A(0) and A* . We calculate that 


Pi = 


Cl — IT , <X> 2 


afd 


Pi =- 


- a—iT~,co 2 


q e 


In addition, according to the equation (3.8), there is 

fO fO — 

< 

J-l Jf =1 


p\s),p{9)> = p\0)p{0)-\ ^ Q p\Z-9)dr,(9)q(Z)dt. 


= D[l + p 2 p 2 * -\jl,pl)0e^ 0 dq(9) 


ypi) 


= D(l + p 2 p 2 + T lm qe ittV2 - PT 2 p 2 p 2 e - iahT ' ), 
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Further, let (q* (s),q(0)} = 1 , there is D = (1 + p 2 p 2 + ^\* c l e u ° lZ2 — Px 2 p 2 p 2 e ,a>2T2 ') 

Next, according to the algorithm in paper [17] and the similar calculation process with the paper [18], following parameters 
expression be getted. 

8 20 = 2 T* 2 D(aq" +p 2 q"e- 2 ”**'-), 
g u = x(b(2aq + p\qe iT ** + e ~ iT ^ ), 
g 02 = 2x 2 D(aq + p 1 qe lTx ' ah ), 

= 2r;b[^(2tC(0) + iw 2 ™(0))-4 A % (2<>(0 )/j 2 +<>(0)ft)], 

M 2 t 2 


where 


W 2O (0) = ^piOy^ + ^^piO)^ 0 + E x e 2ir '^ e , 


t 2 cd 2 


W n (9) = -- 'M^p(0)e tr ^ +^p(0)e-' T ^ +E 2 . 


3i r 2 co 2 


i t\(o\ 9 . l 8u 


-ir^co-yO 


x 2 co 2 


x 2 co 2 


While, E x and E 2 be determined by the following equation. 


(3.9) 


E x = 2 | 


aft 

-qe** 1 Tl 2ico* 2 +S + j3e 


2 io) 2 - a 

—2 iaCr- 


-HO-yT-y 


y 1 

( " 3 

aq 

/ 

y qe~ 2i6il \ 


E 2 = -\ 


a -a /3 

U ~(a + j3) 


2 aq 


IE\*(0\ , l ?i#Oi 


{q e " ' +e 


In the end, we can calculate the following values of the coefficient. 

C 1 (0) = T -^(*2o*n-2|«, l |-i|«„ 3 | 2 ) + 1 f, 

2 x 2 a> 2 3 2 

Re {(2,(0)} 

/2 ReU'(T*)}’ 

Pi = 2Re{C 1 (0)}, 

T Im{C 1 (0)} + // 2 (ImU'(r;)}) 

2 * * 
x 2 (o 2 

Theorem 3.1 

(i) The direction of the Hopf bifurcation is decided by the symbol of p 2 . If p 0 > 0 , then it's supercritical. If p 2 < 0 , it is 
subcritical. 

(ii) If J3 2 >0 (p 2 < 0 ), then the bifurcating periodic solutions is stable (unstable). 

(iii) If T 2 > 0 (T, < 0 ), then the period of the bifurcating periodic solutions increase (decrease). 
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IV. Numerical Simulation 

In this section, we shall conduct some numerical simulations to verify the main conclusions of this paper, and discuss the 
economic meanings of the different dynamic behaviors. In numerical simulations, according to literature [19], we choose the 

parameter. Let y = 0.5, OC = 2.21 5, q =0.4, q = — 0.088, [5 = 0.8 , 5 = 0. 1, then Eq. (2.2) can be written as 


u l = —0.2215u l -1.82 u 2 -0.2202 (Uj) 2 , 

< 

u 2 =0.4z/ 1 (f-r 1 )-0.8w 2 (f-r 2 )-0.088(M 1 (/-r 1 )) 2 -0.1w 2 . 
With calculation by Mathematica software, when Z { = Z 2 = 0 , we get A w =1.3 1992, 


(4.1) 


A n = 1 . 36142 . So the condition ( H l0 ) is satisfied. WhenTj > 0 , 7 "-, = 0 , we get O) l0 = 0.62444 Z" 1{) = 2 . 10366 . If 

e[0,r 10 ), the equilibrium E* of system (2.2) is asymptotically stable. Once Z { more than critical value 
r 10 = 2 . 10366 , the system (2.2) produce a bunch of periodic solutions nearby the equilibrium point E . Numerical 
simulation results are shown in Figs. 1 and 2. Similarly, when Z", = 0 , Z-, > 0 , we get (O 20 = 1 .4739 1, Z 20 = 1.21 49 (Figs. 3 
and 4). When Z"j = Z" 2 = 0 , we get CO 0 = 1 . 1 152 , T 0 = 0.95603 (see Figs. 5 and 6). When T x > 0 , X 2 = 1.2 e ( 0 , r 20 ) , we 
get 0)* = 1 . 3327 , 7 * = 0.2342 (see Figs. 7 and 8). Similarly, when r 2 > 0 , Z x = 1.8 e ( 0 , r 10 ) , we get 
CO * = 0 . 775 , Z* = 0.6286 (see Figs. 9 and 10). 


In addition, fU 2 = 6. 1344> 0, yS, = —121. 1632< 0. , T 2 = 1 1 .2663> 0 By the theorem (3.1), when 
Z", > 0, Z x e (0, r 10 ) , the Hopf bifurcation produced by the system (2.2) is supercritical, and the period of the bifurcating 
periodic solutions increase. 




Figure 2:The equilibrium E becomes unstable and a Hopf bifurcation occurs when z, =2.108> z w 

"2 ul u2 



Figure 3:The equilibrium E is locally asymptotic al ly stable with z 2 = 1 .3 < r 20 
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Figure 4:The equilibrium E becomes unstable and a Hopf bifurcation occurs when r 2 = 1 .3 8 1 > r 
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Figure 6: The equilibrium E* becomes unstable and a Hopf bifurcation occurs when t, = r 2 = 0.9565>r l( 
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Figure 7:The equilibrium E* is locally asymptotically stable with r, = 0. 1 < zv, r 2 = 1.2 
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Figure 8:The equilibrium E becomes unstable and a Hopf bifurcation occurs when 

r, =0.235> r lt ,r 2 =1.2 
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Figure 9:The equilibrium E 1S locally asymptotically stable with 


t 1 = 1.8, r 2 = 0.2 < r 2 
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V. Conclusion 

Since the anticipation of capital stock and its future value are directly interrelated, the government should consider the 
expectation of capital stock in investment decisions at present stage. At the same time, the implementation of past investment 
decisions also need a pregnancy period which leads to production delay. 

In this paper, the main contribution be lied in the following content: the first, we improved the traditional Kaldor-Kalecki 
model with two delays in the gross product and the capital stock. We set up the Kaldor-Kalecki model of differential 
equation with the two delays; The second, we study the stability and Hopf bifurcation. The results indicate that both capital 
stock and investment lag are the certain factors leading to the occurrence of cyclical fluctuations in the macroeconomic 
system. Moreover, the level of economic fluctuation can be dampened to some extent if investment decisions are made by the 
reasonable short-term forecast on capital stock. And finally conduct numerical simulations to prove the conclusions. 

The above arguments are well prepared to the further research work, but there are many undeveloped theory that need further 
explore. The results of this paper can be used as qualitative analysis tool of mathematical economics and business 
administration. 
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